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1. Introduction 

The gradient flow renormalized coupling is a useful observable both in scale setting and in step 
scaling function studies [1,2]. It is easy to determine with high accuracy, however it can suffer from 
large discretization effects [3-5]. While it is possible to systematically remove the lattice artifacts 
by improving the lattice action, lattice operator, and the gradient flow force at the same time, such 
improvement is not always practical [6]. Tree level perturbative improvement developed in ref. [7] 
is easy to implement and appears to be effective at weak couplings. Recently we introduced the “t- 
shift" improvement, an empirical approach to reduce, even remove the f/(a 2 ) corrections of the step 
scaling function, “t-shift" improvement is a non-perturbative method that can be used and appears 
to be effective even at strong gauge couplings. In this talk I review this method and summarize 
results obtained for the discrete /J function in the Nj = 12 and 8 flavor SU(3) systems [8, 9]. I also 
suggest a generalization of the “t-shift" improvement that can be used in scale setting and illustrate 
it using Symanzik flow data on 2+1+1 flavor HISQ configurations [10, 11]. 

In section 2 I give a brief description of the t-shift improved gradient flow coupling, following 
closely Refs. [8, 9]. Section 3 illustrates the t-shift optimization of the ?o gradient flow scale, 
followed by the step scaling function study of the Nf = 12 and 8 flavor systems in sections 4 and 5. 

2. The t-shifted gradient flow running coupling 

The gradient flow renormalized coupling at energy ft = \/St is defined as 

g 2 GF (Li = l/VSt) = -^{t 2 E(t)), (2.1) 

where t is the flow time, E denotes the energy density, and the normalization constant is cho¬ 
sen such that £g F (jU) agrees with the continuum MS coupling at tree level. The t-shift improved 
gradient flow coupling introduced in ref. [8] replaces gQ F (iu;a) with the t-shifted coupling 

SGF(/^’ a ) = —y + Tofl 2 )) , (2.2) 

where To <C t /a 2 is a small shift in the flow time. In the continuum limit To a 1 —> 0 and y Fi| .-(/./) = 
gQ F (/i). At finite lattice spacing gj ]]: (fL) differs from gQ F (fl) by an A(to« 2 ) term. Since gQ F (fl) 
itself has G[cr) cut-off corrections, an appropriate choice of To can remove those, making the 
t-shifted coupling improved 

sipt^a) =gGF(ir,a = 0) + iff(a 4 [loga) n ,a 4 ). (2.3) 

The t-shift improvement contains only a single parameter, To- To achieve full (A (a 2 ) improve¬ 
ment this To must depend on other parameters, like gQ F (fl) and the bare gauge coupling. Fortu¬ 
nately in practice it is sufficient to choose To to be a constant or only weakly gj }]: (fi ) dependent to 
remove most ff{a 2 ) lattice artifacts. The t-shift does not have to be fully optimized, even if To does 
not remove all (A(a 2 ) corrections, the t-shifted coupling has the correct continuum limit. 

The t-shift improved t a lattice scale can be defined as t 2 (E(t a + To)) = oc . The traditional 
to scale correspond to a = 0.3 but any other (reasonable) value can be used. Requiring that the 
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Figure 1 : Left panel: The gradient flow running coupling gQ F { ll ) without t-shift improvement (To = 0 ) of 
the m s /mi = 27 (m3) data set as the function of t jUp Errors are smaller than the width of the lines. The 
dashed lines signal where \[%1 /a < 2 and lattice artifacts could be large. Right panel: Same as the left panel 
but with near-optimal To = — 0.18 t-shift improvement. 


relative scales t a /to are independent of the lattice spacing, i.e. identical on different configuration 
ensembles, defines the optimal To value as will be illustrated in section 3. 

In step scaling function studies the optimal To removes the f7(a 2 ) terms of the discrete (5 
function 


P hll (gf,s-a) 


Sc{L\a) —g^(sL',a) 

log 0 2 ) 


(2.4) 


The subscript c in g~ relates the energy scale and the lattice size, y/St = l/fl = cL [12], its value 
is usually c = 0.2 — 0.4. 

Since the gradient flow is evaluated through numerical integration, the replacement g^ F —>• g^ F 
can be done by a simple shift of t without incurring any additional computational cost both in scale 
setting and step scaling studies. 


3. The t-shift improved lattice scales 

I illustrate the t-shift improved scale setting with 2+1+1 flavor HISQ ensembles [13, 14] using 
Symanzik flow data generated for ref. [10]. I consider the ensembles with physical strange mass 
and three different light to strange quark mass ratios: m s /mi = 5 (ml), 10 (m2) and 27 (m3). The 
ensembles, according to their lattice spacings, fall into four groups, a ~ 0.15fm, 0.12fm, 0.09fm 
and 0.06fm (ensembles A, B, C and D). Further details can be found in ref. [10]. 

The lattice artifacts of to /a 2 for the m s /mi = 27 data set are illustrated in the left panel of 
figure 1 where I compare the gradient flow running coupling gQ F (fj;a) as the function of the flow 
time t normalized by to for the four different lattice spacings ensembles. Since the different curves 
are rescaled with their corresponding to / a 2 lattice scale, they are forced to cross at t/to = 1 . 

The condition t 2 (E) = 0.3 to set the scale is arbitrary, any (reasonable) value of a in the 
relation t 2 (E(t a + To)) = a should predict consistent lattice scales. If the gradient flow scale had 
no lattice spacing corrections, the curves corresponding to different ensembles in the left panel 
of figure 1 would overlap in a wide range, limited only by finite volume effects at large t and by 
gradient flow integration effects at small t. This is obviously not the case. 
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Figure 2: \JtoJt\ as the function of a 2 /to. The different colors/symbols correspond to different m//m s ra¬ 
tios. Left: To = 0, no t-shift improvement. The dashed lines are quadratic extrapolations. Right: To = —0.18 
improved gradient flow. The dashed lines are linear extrapolations excluding the coarsest (A) configuration 
set. 


To quantify the lattice artifacts I compare two different scales, to at a = 0.3 and t\ at a = 0.35. 
The left panel of figure 2 shows the ratio y / to/t\ as the function of the lattice spacing a 2 for all 
three m s /mi data sets. This ratio would be constant if there were no lattice artifacts, but as the 
figure shows \ftoJh is not only not constant, it does not even show a linear dependence on a 2 . 

Introducing a small To shift changes the predicted scales slightly and improves the scaling 
considerably. The right panel of figure 2 shows the ratio \Jto/t\ with To = —0.18. This value was 
obtained by trial-end-error without attempting to identify a precise optimal value: any To is correct 
in the sense that it should predict the same continuum limit. 

Since I use the same To value at every bare gauge coupling, I do not expect full f)‘(cr) improve¬ 
ment. Nevertheless the three finest data sets predict a nearly constant value for sjtojh suggesting 
only minimal lattice artifacts. At the same time the coarsest “A" data set at a k, 0.15fm shows 
significantly cut-off corrections and cannot be included in any reasonable continuum limit extrap¬ 
olation. It appears that the coarse a ~ 0.15fm configurations arc not in the fl/cr) scaling region of 
the quantity to- Without t-shift improvement the rather large lattice artifacts mask this problem. 

The right panel of figure 1 is the analogue of the left panel showing the gradient flow run¬ 
ning coupling now with near-optimal To = —0.18 t-shift for the m3 ensembles. g 2 . ]: (fi;a) shows 
little variation with the lattice spacing after the initial gradient flow integration artifacts die out, a 
significant improvement over the To = 0 case in the left panel. 

It is worth considering other lattice scales to check lattice artifacts with and without the t-shift 
improvement. As an example in figure 3 I show y/to/wo with To = 0.0 and -0.18. Lattice corrections 
are significant without t-shift, while the t-shifted ratio is almost constant and can be extrapolated 
linearly if only the three finer data sets are considered. As before, the coarsest “A" configuration 
sets do not appear to be in the fj/a 1 ) scaling regime. Comparing to to the more conventional r\ 
Sommer scale leads to very similar conclusions. 

In a recent study of 4 light and 8 heavy flavors we used the gradient flow coupling in a similar 
way. We found that a small t-shift removed most observable lattice artifacts in the infrared, showing 
the emergence of walking behavior as the mass of the 8 heavy flavors decreased [15] 


4 


















Improved gradient flow for step scaling function and scale setting 


Anna Hasenfratz 


s 


1.1 

1.05 

1 

0.95 

0.9 

0.85 

0.8 

0 0.2 0.4 0.6 0.8 1 1.2 


T 0 = 0.0 


ml i—B— 

m2 - 

m3 i— *- 






.S''- 

3r 


..S' 
' M 


£ 


1.1 

1.05 

1 

0.95 

0.9 

0.85 

0.8 

0 0.2 0.4 0.6 0.8 1 1.2 


T 0 = -0.18 


. 


—7®—- 

■-~M . 


ml i—B- 

m2 i-= 

m3 i— 


a 2 /t 0 


a 2 1 to 


Figure 3: Same as figure 2 but for the ratio v /7q/vvq. 



Figure 4: Left paneLThe A/- = 12 mnning coupling gfiL) versus the bare coupling /j/ on several volumes, 
for c = 0.2. Right panehContinuum extrapolations of the 12-flavor finite volume IRFP gl(L), with several 
different f-shift coefficients To for fixed scale change s = 2. 


4 . Nf— 12 flavors 

Thel2-flavor 3-color system has been investigated by several groups with somewhat contro¬ 
versial conclusions. While recent finite size scaling analysis suggests that the infrared behavior of 
this system is consistent with conformality and mass anomalous dimension y„, = 0.235(5) [16, 17], 
a step scaling function study remains the most reliable method to identify its infrared fixed point. 

In ref. [8] we used nHYP smeared staggered fermions and considered six different volumes up 
to 36 4 . We performed the simulations in the m = 0 chiral limit and the range of volumes allowed 
us to carry out step scaling analyses with scale changes 5 = 4/3, 3/2 and 2. We chose the gradient 
flow coupling with c = 0.2 to minimize the statistical errors, though we verified that other c 
values (0.25 - 0.35) gave consistent results. 

The left panel of figure 4 shows the running coupling gf (L) as the function of the bare gauge 
coupling ffl- for different volumes. The interpolating curves are polynomial fits that cross in the 
range 6.0 < Pf < 6.5. The crossing from lattices with linear size L and sL defines the finite- 
volume IRFP coupling g+(L;s). If the IRFP exists in the continuum limit then the extrapolation 
1 im /0 gf(L;s) = gl has to be finite and independent of the scale change s. The right panel of 
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Figure 5: Continuum-extrapolated discrete fl function for scale change s = 3/2 with c = 0.25 (left) and 
0.3 (right). In each plot we include once- and twice-smeared results using the optimal T op t = 0.07 and 0.18, 
respectively, as well as two-loop perturbation theory (solid line) and the four-loop perturbative prediction in 
the MS scheme (dotted line). The darker error bands indicate our statistical uncertainties, while the lighter 
error bands show the total uncertainties, with statistical and systematic errors added in quadrature. 


figure 4 illustrates the continuum extrapolation of g1(L) with scale change s = 2 for various choices 
of the (-shift parameter To. Most lattice cut-off effects are removed with (-shift, T op t ~ 0.04. 

The near-optimal T opt ~ 0.04 turns out to be also near-optimal for scale changes s = 3/2 and 
4/3, making the extrapolation to the continuum very stable, predicting that the IR fixed point is at 
renormalized coupling g 2 = 6.18(20) in the c = 0.2 scheme. 


5. Nf = 8 fundamental flavors 

The Nf = 8 flavor SU(3) system has received significant attention as a chirally broken system 
near the conformal window and candidate model for walking behavior. 2-loop perturbation theory 
predicts that Nf = 8 is below the conformal window though at 3- and 4-loop the MS /3 function 
develops an infrared fixed point at strong coupling. In ref. [9] we studied the step scaling function 
of this model using two different lattice actions, one with once the other with twice nHYP smeared 
staggered fermions. The continuum extrapolated step scaling function must be independent of the 
lattice action so comparing different actions at identical renormalized couplings can give insight of 
systematical errors. 

Our numerical simulations were carried out in the m = 0 chiral limit. We did not observe 
spontaneous chiral symmetry breaking even though we probed fairly large renormalized couplings, 
reaching g%(L) = 18.0(1) on our largest 30 4 volume with two nHYP smearing steps. 

Figure 5 summarizes the results for the continuum-extrapolated discrete /3 function both with 
c = 0.25 and 0.3. The once- and twice-smeared results are consistent and predict a /3 function that 
is significantly smaller than the perturbative universal 2-loop value. Our results are even below the 
(non-universal) 4-loop MS prediction, though the non-perturbative results increase monotonically 
in the range accessible in our simulations. 

It would be very interesting to repeat this calculation with Wilson or domain wall fermions. 
Those actions might be able to probe stronger couplings and reach the expected chirally broken 
regime. Testing universality between different fermion formulations would be important as well. 
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6. Conclusion 

In this paper I described several applications of the “t-shift" improved gradient flow coupling. 
I illustrated this approach for scale setting on 2+1+1 flavor HISQ configurations and showed that 
after optimization the ratios ^/ip/wo or y/to/ri show negligible dependence on the lattice spacing. 

The t-shift improved coupling reduces ff(a 2 ) lattice artifacts of the discrete /J function, im¬ 
proving the accuracy of continuum extrapolations in step scaling studies. Our results in the Nf = 12 
SU(3) system predict an infrared conformal fixed point. With 8 flavors we could follow the discrete 
/3 function up to g 2 ( L) ~ 14. We found that the discrete /3 function, while significantly smaller than 
the 2-loop prediction, increases monotonically up to this gauge coupling. 
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